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Abstract 

Potential (electrostatic) surface waves in plasma half-space with degenerate electrons are studied 
using the quasi-classical mean-field kinetic model. The wave spectrum and the collisionless damping 
rate are obtained numerically for a wide range of wavelengths. In the limit of long wavelengths, 
the wave frequency u approaches the cold-plasma limit u = u p /y/2 with oj p being the plasma 
frequency, while at short wavelengths, the wave spectrum asymptotically approaches the spectrum 
of zero-sound mode propagating along the boundary. It is shown that the surface waves in this 
system remain weakly damped at all wavelengths (in contrast to strongly damped surface waves in 
Maxwellian electron plasmas), and the damping rate nonmonotonically depends on the wavelength, 
with the maximum (yet small) damping occuring for surface waves with wavelength of ~ 57rAi?, 
where is the Thomas- Fermi length. 
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I. INTRODUCTION 



It has long been known [l-3] that bounded plasmas support a special type of collective 
electrostatic and electromagnetic excitations - the surface plasma waves - whose field and 
energy are concentrated near, and propagate along, the plasma boundaries. The surface 
waves (SW) in various (classical) plasmas have been extensively studied, and they have 
found many applications (see Ref. |4| and references therein). Yet recently, due to a re- 
markable progress in nanotechnology, the interest in surface waves supported by various 
nanostructures (especially metallic structures such as thin films and tiny metallic particles), 
and in their interaction with light, has been revived. It is believed that light-induced surface 
excitations in such structures may offer a route to faster, smaller, and more efficient electron- 
,e 8 , a 8 wen as new te C n„o.o gy ff. In p arti cu,„ one cou.d note suen recent advents in the 
new and promising area of quantum nanoplasmonics as the development of the concept of 
spaser Q followed by its further development into a lasing spaser 7], and the experimental 
demonstration of a spaser-based nanolaser js, 9]. 

In view of these developments, understanding the properties of surface waves in various 
metallic (and semiconductor) structures, bounded by vacuum or dielectric, is thus impor- 
tant. Such understanding requires using models describing the dynamic response of such 
structures to self-consistent electromagnetic fields, that also appropriately take into account 
the relevant quantum effects arising from quantum nature of free charge carriers in such 
structures and, in general, from their (quantum) interaction with each other and with the 
underlying ion lattice. These quantum effects may significantly alter the properties of the 



surface waves; see, e.g., Refs. 
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Recently, the dispersion relation of surface waves in one of the basic structures modeling 
a nanoplasmonic device - a semi-bounded collisionless quantum plasma with degenerate 



electrons - was obtained in Ref. 



121 ] using the quantum fluid theory (QFT) approach 13 ]. 



In the electrostatic limit, the authors of Ref. [12[ obtained for the frequency of surface waves 
in this structure: 

!i =^( 1+ vf jr, \A T ^f)- w 

where f2 = co/uj p , Ku = UuXf, H = hu) p /2mvp, co p = (A r Ke 2 n e /m e ) l l 2 is the electron plasma 
frequency, vp = h\/3ir 2 n e /m e is the electron Fermi velocity, = vp/y/3u p is the Thomas- 
Fermi length, to is the SW frequency, k\\ is the SW wave vector along the plasma boundary, e, 



m e and n e are electron charge, mass and number density, respectively, and h is the reduced 
Planck constant. However, the validity of the dispersion relation ([1]) obtained in Ref. [l^J 
(as well as of the similar dispersion relation obtained in Ref. 11]) is limited by the validity 
of the QFT approach itself and is restricted only to long waves, K\\ <C 1. Moreover, 

the QFT approach by its nature completely ignores the purely kinetic effect of collisionless 
damping of surface waves, which is known to be significant, e.g., for potential SW in warm 



electron plasma at short wavelengths 15]. To overcome these limitations, a kinetic approach 
to the problem of SW in semi-bounded plasma with degenerate electrons is needed. 

In this paper, we study potential surface waves in semi-bounded collisionless quantum 
plasma with Fermi-degenerate electrons, using the initial value problem solution for the semi- 
classical Vlasov-Poisson system Q . We obtain the dispersion and collisionless damping rate 
of these waves, valid for both long and short wavelengths, and report on a surprising result 
that these waves remain weakly damped for all values of Ku (i.e., for all wavelengths), unlike, 
e.g., in plasma with Maxwellian electrons in which the SWs are weakly damped only for 
K\\ <C 1 and quickly become strongly damped as K\\ increases. We also report a nontrivial 
nonmonotonic dependence of the damping rate on Kn, featuring a maximum damping at 
K\\ 0.4 (corresponding to the surface wave length of ~ birXp). 

II. METHOD 

A. Model and Assumptions 

We consider a uniform nonrelativistic quantum plasma consisting of mobile electrons 
with charge e, mass m e and number density at equilibrium n e = Uq, and immobile uniform 
background of singly charged ions with number density Uq that neutralizes the electron 
charge at equilibrium. The plasma is assumed to be confined to the region x < by a sharp 
perfectly reflecting wall located at x = 0. We will be interested in evolution of an initial 
perturbation to the equilibrium state of the system's electronic component applied at time 
t = 0. 

In general, this system of many interacting quantum particles (electrons) can be de- 
scribed by the quantum analog of the Bogoliubov-Born-Green-Kirkwood-Yvon (BBGKY) 



hierarchy 
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17] of equations for the electron j-particle quantum distribution functions 
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Wj(Tj,pj,i) (also called the j-particle Wigner functions), where j = 1,2, ...,N, Tj and pj 
are the 3j'-dimensional vectors denoting the sets of coordinates and canonical momenta of 
system particles, and N is the total number of electrons in the system. In this hierar- 
chy, each of the equations for the j-particle quantum distribution function contains the 
(j + l)-particle distribution function, making the whole set of N equations coupled, and 
thus prohibitively large to solve. In practice, however, this hierarchy can be truncated 
by making a physically justified assumption about correlation of particles due to their in- 
teraction. In particular, for a system of weakly interacting particles, with a small plasma 
coupling parameter Q = U int /e khl <C 1 (here U int is the characteristic potential energy of par- 
ticle interaction, and e k[n is the characteristic kinetic energy of plasma particles), the two- 
and higher-order particle correlations can be ignored, leading to the collisionless mean-field 



approximation (18| involving only one equation for the one-particle quantum distribution 
function Wi(r, p, t), where r and p are now the 3-dimensional vectors of particle coordinate 
and momentum. In the quasi-classical approximation, with the effect of quantum recoil 



ignored (see Sec. IIII CI) , this equation reduces to the Vlasov equation for the one-particle 



classical distribution function /(r, v, t), where v = p/m e is electron velocity. For a system 
of particles with electrostatic interaction, in the chosen geometry, the Vlasov equation for 
/(r, v, t) = f(x,r\\,v x ,v\\,t) (where rn and vii are, respectively, the components of r and v 
parallel to the boundary) reads 

dt x dx " dr\\ m e \dx dv x dr\\ dv\\) 
where the electrostatic potential 4>(x,T»,t) is defined by the Poisson's equation 



Aire 



J /(r, v,t)rf 3 v - n 



(3) 



In the absence of fields, the equilibrium distribution function of plasma electrons /o( v ) is 
defined by the Pauli's exclusion principle for fermions, and for low electron temperatures 
Te/^F 1 (where T e is the electron temperature, ep = m e v F /2 = (37r 2 ft 3 n e ) 2//3 /2m e is the 
electron Fermi energy) it becomes 

ISJ _ (Mp (4) 
if v > Vp, 

corresponding to fully degenerate electron distribution. 



The condition of specular reflection of plasma electrons off the boundary at x = implies 
f(x = 0,r||, -v x ,\\\,t) = f(x = 0,r\\,v x , Vy,t). (5) 

B. Initial value problem 

We now introduce a small initial perturbation f p (x, ry, v x , Vy, t = 0) to the equilibrium 
electron distribution function fo(v), \ f p (x, r\\,v x , vy, t = 0)| <C fo( v ), an d study the resulting 
evolution of the system's charge density p(x, ry , t) — e [J f(x, ry , v, t)d 3 v — no] , and hence of 
the electrostatic potential <fr(x, ry,t) defined by Introducing the dimensionless variables 
f2 = to/ Up, K = kAi?, V = v/vf, X = x/ Xf, Ry = ry/Ai?, A^ = vf/VSu p , and following 
Guernsey jsj, the solution of the formulated initial value problem for p(X, Ry,T) with the 
boundary condition (JSJ) is 

p(X,Ry,T) = en p(X,Ry,T), (6) 

where 



-oo 



p(X, Ry, T) = — -- y dff, y rf 2 Ky e 4K n R ii p k (T), (7) 

i pia -\-oc 

p K (T) = — / p(fi,K) e" inT dfi, with a > 0. (8) 

2 71 " J io-oo 

The integration in (jSJ) is performed in complex Q plane along the horizontal contour that 
lies in the upper half-plane Im(fi) = a > above all singularities of the function p(f2,K). 
The function p(f2,K), defined as the Laplace transform of pk(T): 

POO 

p(fi,K)= / p K (T) e iQT dT, (9) 
Jo 

is given by 

G(V,K) 



,5(0, K) = — — / cfV 



e(n,K)J n-V3K-V 



iK\\ 

+ 



1- 1 



^3 V g(y,KQ , (10) 



K l2 e{Q,K')J fi- v / 3K'-V ; 



where the Fourier transforms G(V, K) and G(V, K') of the (dimensionless) initial pertur- 
bation are defined by 

/+oo 
dXe^^K.Vi.Ki), (11) 
-oo 

/+oo 
dX e~ iK * x g(X,V x ,V h K {l ), (12) 
-oo 



with 

g{X,V x ,V h K ]{ ) = J d 2 R,| e" iK "- R " f p (X, R„, V x , V„, 0), (13) 

/ P (X,R||,T4,V||,0) = ^/ P (X,R||,T4,V||,0), (14) 

n 

where 7f|| = |K|||, K = |K|, K = (K X ,K\\), K' = |K'|, and K' = (7C,K||). The functions 
e(Q,K) and C(0,7C||) in (HHD are defined (for Im(O) > 0) as follows: 

1 [K-df (V)/dV^ 



with 



e(Q,K) = 1 - -= — / cfV, (15) 

V ' V3K 2 J Q-V3K-V V 7 

1 K\\ f +co dK 



/o(V) = ^/o(V) = ^/ (v)U fV . 



n n 

For fully degenerate plasma with electron distribution (jlj), the function e(0, K) becomes [15 . 
1SJ: 



o , /o + v^tA 

1 — — In 1 



(17) 



2V3K 11 ^O- / 

where ln(z) is the principal branch of the complex natural logarithm function. 

Note that the solution (fTUI) differs from the corresponding solution of the transformed 
Vlasov-Poisson system for infinite (unbounded) uniform plasma only in the second term 
involving C(0, K\\); indeed, this term appears due to the boundary at x — 0. 

The definition ([9]) of the function p(0, K) of complex O has a sense (i.e., the integral in 
([9]) converges) only for Im(O) > 0. Yet the long-time evolution of p\a(T) is obtained from 
(JHJ) by displacing the contour of integration in complex O plane from the upper half-plane 
Im(O) > into the lower half-plane Im(O) < [^(J. This requires the definition of p(0, K) 
to be extended to the lower half-plane, Im(O) < 0, by analytic continuation of (fTOj) from 
Im(O) > to Im(O) < 0. Hence, the functions 

7(0, K)= [d'V G{V 1 K) , (18) 

e(Q,K), and £(0,7^) that make up the function p(0,K), must also be analytically con- 
tinued into the lower half-plane of complex O, thus extending their definition to the whole 
complex O plane. With thus continued functions, the contributions to the inverse Laplace 
transform (JSJ) are of three sources jsj: 
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1. Contributions from the singularities of J(f2,K) in the lower half of complex Q plane 
(defined solely by the initial perturbation G(V,K)); with some simplifying assump- 
tions about the initial perturbation js| these contributions are damped in a few plasma 
periods and can be ignored. 

2. Contribution of singularities of l/e(Q,K) in the lower half of complex Q plane, of 
two types: (i) residues at the poles of l/e(Q, K), which give the volume plasma os- 
cillations jsj], and (ii) integrals along branch cuts (if any) of l/e(Q,K) in the lower 
half-plane of complex Q, which can lead to non-exponential attenuation of the volume 
plasma oscillations 



21 



22j. 



Contribution of singularities of K\\) in the lower half of complex Q plane, of 

two types: (i) residues at the poles of K\\), correspondingly the surface wave 

solutions of the initial value problem in the considered system 3j, and (ii) integrals 
along branch cuts (if any) of K\\) in the lower half-plane of complex Q, which 

will be discussed elsewhere. 



In this paper, we will only consider the surface wave solutions due to the contribution of the 
residues at the poles of 1/C(^, K\\)- The dispersion and damping properties of these surface 
wave solutions are defined by the dispersion relation for plasma surface waves j^] 

m,K l{ ) = 0, (19) 

which in case of stable plasma (with velocity distribution of electrons having only one max- 



imum) only has non-growing solutions Q = Q{K\\) G C with Im(fi) < 23] 



C. Weakly damped surface waves 

Out of all complex- valued solutions of f[T9l . Q(K\\) = Q s (Kn) + iT s [K») (here Q s € R is 
the dimensionless frequency, and T s < is the dimensionless damping rate of the surface 
wave), only those corresponding to weakly damped surface waves, with |r s (iC||)/f] s (ifii)| 1, 
are of physical interest. Such solutions can be obtained by solving, instead of the actual 
dispersion equation (fl9~l) . the following set of approximate equations (that follow from (fl9~|) 
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for |r s (^||)/a(^||)|<l): 



Re [C(n a ,K\ 



0, yielding f2 s = Q S (K\\) e 
Im [C(fi s ,^||)] 



d Re [C(fi a ,AT||)] /dVt s 



(20) 
(21) 



n s =n s (K||) 

which only involve £(f2, ify) as a function of real Q. Thus, instead of performing analytic 
continuation of £(Q,K\\) from Im(f2) > to Im(fi) < 0, it suffices to only analytically 
continue £(Q,Kn) from Im(Q) > onto the real axis Im(fi) = of the complex Q plane, 
thus defining it for AgR. This analytic continuation of ((fl, K\\), required for the numerical 
solution of Eqs (I20]) - (I2~T]) . is discussed in Appendix lAl 



III. RESULTS AND DISCUSSION 

The numerical solution of Eqs ( 120]) and f l2~T]) for a given K\\ > 0, with ((Q,K\\) defined 
for Q G K as discussed in Appendix |Al yields the dispersion Q S (K\\) and damping T 3 (K\\) of 
surface waves; below we discuss them in greater detail. 



A. Dispersion 

The dispersion curve of surface waves Q S (K\\) is shown in Fig. [Q The SW frequency 
increases monotonically with Kh, starting from the well-known long-wavelength limit of 
fl s = l/v2 at K\\ — > (which also corresponds to the cold-plasma limit of SW frequency), 
and increasing as 



n s (K\\ 



i 



(1 + 0.95^1) 



(22) 



for K\\ 1. We note the discrepancy between (122]) and the corresponding small- K 
asymptote of Lazar et al.'s result (JTJ of Ref. [12|: Vt s {K\\) w (l/v^)(l 



Z/2K\\ 



[l/y/2){l + 1.23^1). This discrepancy can not be attributed to the effect of quantum 



recoil, which is ignored here but included in the QFT model of Ref. [12], as the quan- 
tum recoil only gives the higher-order contribution oc H 2 K^ to the dispersion (pQ) at small 
K\\. In fact, the mentioned discrepancy is due to the error in the coefficient of the clas- 
sical pressure gradient term used in Ref. 12j: instead of the three-dimensional pressure 
gradient term (3/5)vpV(n e /rio) 24|, |25|, the one- dimensional term VpV(n e /rio) was used. 
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FIG. 1: (Color online) The dispersion curve of surface waves £l s (K\\) (solid blue line), and the 
asymptotes: the long- wave limit l/\/2 (dotted line), the long- wave asymptote ([22]) (dot-dashed 
blue line), and the short-wave asymptote ([23]) corresponding to zero sound propagating along the 
boundary (dashed red line). 



With this error corrected, the QFT model of Ref. 12] yields the small-i^y asymptote of 
Q S (K\\) = (1/V2)(1 + y/9/10K\\) that matches our asymptote ([22]) . 

At short wavelengths, K\\ 3> 1, the SW dispersion can be approximated as (see Ap- 
pendix [B]) 

VL S « y/$K\\ (1 + 2 exp [-2 - 4fif|] ) , (23) 
and quickly approaches the line VL = V3 K \\ t hat corresponds to the dispersion of volume zero 



sound mode in degenerate Fermi gas 19 



26j propagating along the boundary, with K = Ku . 
Note that for K\\ ^> 1 the SW frequency Q s becomes practically indistinguishable from the 
corresponding frequency Q v of the volume plasma wave propagating along the boundary 



(with K = Ku 



1 2 1 



V3K\\ (l + 2 exp [-2 - 2Kf\) . 
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B. Damping 



The K\ | -dependence of the surface wave damping rate r a (iT||) is shown in Fig. [2j In the 
long-wave limit K\\ <C 1, the damping rate increases linearly with K\\ as 

|r s (A'n)| « 2.1^3 • 10~ 2 K h (24) 

similar to the damping rate of SW in plasma with Maxwellian electrons j^, [lsj]. However, 
as K\\ increases, an important difference between SW damping in plasma with degenerate 
electrons and in plasma with Maxwellian electrons becomes obvious. In Maxwellian plasma, 
the SW damping rate increases monotonically with Ku , quickly exceeding the SW frequency, 
so that the SW become strongly damped at short wavelengths {3, 15]. In degenerate electron 
plasma, however, the SW damping rate has a non-monotonic dependence on K\\, as seen in 
Fig. [2J at small K» , the damping rate increases almost linearly with K» , reaching a distinct 
maximum of \T S \ « 6.2 • 10~ 3 at K\\ ~ 0.4, and then decreases monotonically with K\\ for 
K\\ > 0.4, approaching zero at large Kn. 

Most importantly, as seen from Fig. [21 the SW damping rate in degenerate electron 
plasma remains small (|r s /f2 s | <C 1) at all values of K\\, and hence the electrostatic surface 
waves in degenerate collisionless electron plasma are weakly damped at all wavelengths (yet 
the surface waves with K\\ = fcnA_F ~ 0.4 are preferentially damped, as compared to other 
wavelengths). 

The following note on the nature of damping of surface waves in collisionless plasma is 
in order. It is well known that in a Maxwellian plasma the collisionless (Landau) damping 
of surface plasma waves is not exponentially small and significantly exceeds that of the 
volume plasma waves in unbounded plasma, even at small value of K\\ when their phase 
velocity greatly exceeds the electron thermal velocity {3} [lsj. The reason is that, as seen 
from Eq. (fT6l) , the surface wave field is essentially a result of plasma response to a collection 
of "virtual" plasma waves [i.e., having frequency Q and wave vector K, but not having 
a dispersion relation Q = Q(K), unlike the "real" plasma waves with dispersion relation 
Q = Q(K) following from e(Q, K) = 0], with the same frequency Q = Q s but with different 
wave vectors K with absolute values K = \J + ranging from K = K\\ (for K x = 0) to 
K — > 00 (for K x — > 00). Each of these virtual waves with its own K interacts resonantly with 
plasma electrons whose velocities along K is close to the phase velocity of the virtual wave 
(and such electrons always exist in Maxwellian plasma), and thus is subject to a non-zero 
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FIG. 2: (Color online) The absolute value of surface wave damping rate |r s (ET||)| (solid blue line), 
and the long- wave asymptote (|2T|) (dashed line), as functions of K\\. 

Landau damping. As a result, the surface wave, consisting of both weakly damped virtual 
waves (with Q < K, where for Maxwellian plasma K = kXp, Xd = VTe/^p is the electron 
Debye length) and strongly damped virtual waves (with Q > K), interacts resonantly with a 
significant part of the plasma electron distribution, and hence is strongly damped, compared 
to the "real" plasma waves with the same wavelength. 

In the degenerate plasma considered here, the difference in damping of surface and the 
"real" volume plasma waves is even more striking. Indeed, the "real" volume waves in such 
plasma, defined by the dispersion relation e(Q,K) = with e(Q,K) given by Eq. (fTTj) . are 
not subject to Landau damping at all, as their phase velocity exceeds the maximum possible 
velocity of degenerate plasma electrons, f2 > \[?>K (here K = k\ F , Xp = v F /y/3u p is the 
electron Fermi length), for all values of K 19[. The "virtual" plasma waves are, however, 
still subject to Landau damping by plasma electrons if their phase velocity is less than the 



electron Fermi velocity, Q s < V3K, i.e., if their K x is large enough, K x > — Kp 

such virtual waves always give a contribution to the surface wave field, as K x in ( 1T61) attains 
arbitrarily large values. This also results in the finite, albeit small (yet infinitely large 
compared to the zero damping of the "real" volume plasma waves) Landau damping of 
surface waves in degenerate plasma. 
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C. On applicability of the used approximations 



The collisionless mean- field approximation used here is justified for a system of weakly in- 
teracting particles, when two-particle correlations (and higher-order particle correlations as 
well) can be ignored. Physically, this corresponds to a system in which the collective effects 
dominate over the effects of particle collisions, which happens when the plasma coupling 

1 /3 

parameter Q = t/ int /e kin is small (see Sec. Ill Aj) . With U- mt ~ e 2 n and e^n ~ e F for degen- 



erate electrons, this implies that the collisionless mean-field approximation for degenerate 
electron plasma is justified when 

Q~^~(^V«1. (25) 
Moreover, when the condition (125]) is satisfied, the effect of quantum recoil on the dispersion 



properties of electrostatic plasma waves is also negligible, at least for K < 1 13|, |27j , and 
the only quantum effects come from the Fermi-Dirac statistics of plasma electrons. Thus 
the mean-field quasi-classical kinetic model (|2]), in which particle correlations and quantum 
recoil are neglected, is internally consistent under the condition (125]) . We note that for 
degenerate electron plasma, the coupling parameter Q ~ e 2 n^ 3 /eF scales as n 1//3 and 
thus decreases as the plasma density no increases; hence the condition (125!) is satisifed for 
sufficiently dense plasma, with n ^> (e 2 m* / h 2 ) 3 , where m* is the effective electron mass in 

1/3 

such plasma. However, as Uq increases, the electron Fermi velocity increases as n and may 
become comparable to the speed of light c, in which case the relativistic effects may become 
important; to avoid this, we also require that n <C (m*c/h) 3 . Thus the nonrelativistic 
mean-field quasi-classical approximation (T5]) is valid for quantum plasmas with degenerate 
electrons with densities 

2*\3 / * \ 3 

e m \ m c\ . 

_)«„„«(__). ( 26) 

The condition (|26l) may be satisfied for some semiconductors, as pointed out in Ref. 22]. 

For electron gas in metals, however, the condition ( 125]) (as well as the first part of the 
condition (1261) ) is not satisfied (in fact in metals Q > 1), which suggests that in metals the 



collisionless mean-field approximation (as well as the QFT model used in Ref. 12], that is 



derived from the mean-field kinetic model 



131 ] ) is formally not justified, as the two- and 



higher-order particle correlations (due to their Coulomb and exchange interactions) become 
important. Instead, a kinetic model of moderately coupled plasma, accounting for particle 
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correlations and quantum recoil, should be developed - a difficult task, given the lack of a 
small parameter characterizing the electron interactions in metals. Nevertheless, there are 
reasons, discussed below, to believe that our idealized model can still work well (at least 
qualitatively) in metals, despite being formally unjustified there. 

In the experimental work of Watanabe [28| , the relation between energy loss and scatter- 
ing angle of 25 keV electrons passing through thin metallic films was measured, and it was 
found, in particular, that the value of energy loss increases with the scattering angle of elec- 
trons. An empirical formula representing this relation for Be, Mg, and Al films was found to 
be in good agreement with the corresponding theoretical formula that follows from the dis- 
persion relation of electron plasma waves (excited by the energetic electrons traversing the 
metallic film), derived using the mean- field kinetic model of free electrons in metal, despite 
this model being formally unjustified in metals due to violation of (125]) there. Moreover, the 
effect of quantum recoil, although detectable, is small compared to the effect of Fermi-Dirac 
statistics in the whole range of measured scattering angles, especially at small scattering an- 
gles (corresponding to small values of plasma wave number). These facts can be perceived 
as an experimental evidence supporting the claim that the mean-field quasi-classical kinetic 
model can work well even in metals, despite there being Q ~ 1. 

There is some theoretical evidence of this as well. In Ref. 29[, the effect of electron 



correlations due to their exchange interaction was taken into account under the assumption 
of characteristic plasma perturbation wavelengths being large compared to the correlation 
length, using the Hartree-Fock approximation. In particular, it was found that the exchange 
correlations affect the spectrum of volume plasma waves by modifying their dispersion rela- 
tion (in the limit of small K), which becomes: 

ft 2 = 1 + -K 2 - —QK 2 , (27) 
5 80 ^ ' V ; 

where the third term on the right is due to the exchange correlations, while the second term 
is due to the quantum statistics resulting in the Fermi-Dirac distribution fl4]) of electron 
velocities in degenerate plasma. We see from ( 12T|) that the term due to exchange correlations 
is small compared to the term due to quantum statistics even for Q ~ 1. This suggests that 
in metals the exchange correlations of electrons only lead to minor modification of plasma 
wave spectra, and neglecting them does not lead to a serious error, while greatly simplifying 
the model. As for the quantum recoil - it modifies the spectrum fT27|) by adding the higher- 
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order term (1/16)QK A . Even for Q ~ 1, the term due to quantum recoil remains small 
compared with the term due to quantum statistics for K < 1, and can be safely neglected. 



Finally, it has beem pointed out in Ref. 13| that in plasmas with low electron temper 



atures, T e /ep <C 1 (e.g., for metals at room temperature), the effect of electron-electron 
Coulomb collisions is negligible. In fact, the typical electron-electron collision frequency v ee 
for metals at room temperature is of order v ee ~ 10 10 s _1 which is many orders of 
magnitude smaller than the typical plasma oscillation frequency u p ~ 10 16 s -1 ; hence the 
electron-electron collisions are not expected to play a significant role for processes occuring 
at the characteristic collective plasma time scale r p ~ u~ l (such as the surface waves stud- 
ied here, with u s ~ uj p ), and can thus be neglected. Indeed, the dimensionless surface wave 
damping rate due to the electron-electron collisions, T e s e = ^ e s e /u p with ^ e s e ~ v ee , is of the 
order T e s e ~ 10 -6 , which is much smaller than the characteristic collisionless damping rate 
T s ~ 10~ 3 obtained in this work (see Fig. [2]), and thus can be safely neglected. 

Beside colliding with each other, the electrons can also collide with the ions of metal lat- 
tice, yet these collisions were also neglected in our model. To justify this approximation, let 
us estimate the electron-ion collision frequency v e i in typical metals used in plasmonic appli- 
cations, and compare the characteristic surface wave damping rate Tf associated with these 
collisions with the collisionless damping rate T s ~ 10~ 3 obtained above. The experimentally 
measured electric resistivity of metals such as gold and aluminium at room temperature is 
of order p ~ 1CT 8 Ohm ■ mjso]]. Using the definition of resistivity p = E/j, where E is the 
applied electric field, and j is the current density in metal (which is assumed to be entirely 
due to the free electrons of conductivity), and the Ohm's law j = {eoui/v e i)E for electrons, 
where Eq = 8.8542 • 10 -12 F m _1 is the dielectric permittivity of vacuum, we have (in SI 
units): 

v ei ~ 5 ^pP. (28) 

For metals, u p ~ 10 16 s _1 , and (|28|) yields v ei ~ 10 13 s -1 , which is small compared to the 
characteristic surface oscillation frequency, and thus is not expected to have a significant 
effect on the dispersion of surface waves. On the other hand, the dimensionless surface wave 
damping rate due to electron-ion collisions, Ff = /cu p with ~ u ei , is of order Tf ~ 
u e i/u p ~ 10 -3 , which is comparable to the collisionless damping rate T s ~ 10~ 3 obtained 
here (recall that the maximum collisionless damping rate at K\\ ~ 0.4 is T s pa 6.2 • 10~ 3 ). 
Thus the electron-ion collisions, neglected in our model, may lead to additional collisional 
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damping of surface waves that is comparable to the collisionless damping obtained from the 
collisionless kinetic model. 

The above discussion suggests that the quasi-classical mean-field kinetic model used here 
is adequate for describing surface plasma waves even in moderately coupled plasmas with 
Q ~ 1 such as the electron plasma in metals. To verify this, as well as to assess the qualitative 
and quantitative importance (or otherwise) of the effects associated with moderate values 
of the coupling parameter Q, a comparison of the model's predictions with experiments in 
various quantum plasmas, including the electron plasma in metals, is needed. 

IV. CONCLUSION 

In this paper, electrostatic surface waves in semi-bounded plasma with degenerate elec- 
trons were studied using the nonrelativistic collisionless mean-field kinetic model. The dis- 
persion relation for the waves is obtained from the initial value problem, and its solution 
corresponding to weakly damped surface waves is presented, yielding dispersion and col- 
lisionless damping of the waves for an arbitrary wave number K\\. It is shown that the 
surface waves in the semi-bounded plasma with degenerate electrons are weakly damped at 
all wavelengths, and their damping rate exhibits nonmonotonic dependence on K», linearly 
increasing with K\\ at K\\ <C 1, then reaching maximum at K\\ ~ 0.4, then falling off rapidly 
to zero as K\\ increases. This is in contrast with the strong damping of surface waves in 
semi-bounded plasma with Maxwellian electrons, and is the consequence of the effect of 
quantum statistics (leading to Fermi-Dirac velocity distribution) for plasma electrons. This 
work, using the more general kinetic model, extends the results of Ref. [3] obtained using 
the quantum fluid theory, in two ways: (i) the range of the dispersion relation Q S (K\\) is 
extended from K\\ <C 1 in Ref. [12I to K < 1 in this paper, and (ii) the collisionless damping, 
absent in the model used in Ref. [3], is obtained and discussed here. 
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Appendix A: Analytic continuation of £(Q,K\\) 



To perform the required analytic continuation of £(fi, Kn) onto the real axis of the complex 
Q plane, we start from the definition fTTB]) of ((Q,K\\), with e(Q,K) given by (|17p . for 
Im(f2) > (where £(Q,K\\) is an analytic function of Q), and then reduce Im(fi) down 
to zero, taking the limit Im(f2) — > 0+ while ensuring that the analyticity of ((Q,K\\) is 
preserved in the process. For Im(f2) > 0, C(^? -^||) is defined in terms of the integral 

dK r 



K 2 e(n,K)' 



(Al) 



where K» > 0, K x G R, K = JK 2 X + Kt, and e(fi,K) is given by Eq. (ED. The function 



[K 2 s(p,,K)\ 1 under the integral in flAlh is an elementary function of K = yK 2 + K 2 , 
which in turn is an elementary function of K x . Thus the function [K 2 e(Q, K)] 1 can be 
extended to complex K x plane by analytic continuation from the real axis \m.(K x ) = of the 
complex K x plane, which is achieved by taking the principal branches of the complex square 



root function K = yK 2 + K 2 and of the complex logarithm function ln[(f2 + y/3K)/ (f2 — 
y/3K)] in (fT7|) . considered as functions of complex K x . The resulting function [K 2 e(Q, K)\ 1 
of complex K x has the following singularities in the complex K x plane: 



I. Branch cut of the complex square root yK 2 + K 2 , taken along the negative real axis of 
the argument K 2 + K 2 . This branch cut maps into two branch cuts of [K 2 e(f2, K)\ 1 
in the complex K x plane, given by two parametric equations: 



K x = ±iJK 2 + r, with K\\ > 0, r G [0, +oo). (A2) 



2. Branch cut of the complex logarithm in ( 1TT|) . taken along the negative real axis of 
the argument (f2 + \/3K)/(il — \f?>K). This branch cut maps into two branch cuts of 
[K 2 e(f2, K)\ 1 in the complex K x plane, given by two parametric equations: 



K x = ±iJ y (7^) ~ K \\ ' with 0eC - K \\ > 0, r G [0, +oo). 



(A3) 



3. Two poles K x = ±iK\\ (K\\ > 0) at the roots of K 2 = 0, lying symmetrically above 
and below the real axis of the complex K x plane. 

4. Two poles ±K r x 6 Cat the roots of K) = 0. Note that for any K G R^effi, K) = 



does not have roots with Im(fi) > 0, if the plasma equilibrium is stable [23|, which is 
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FIG. 3: (Color online) Singularities (branch cuts and poles) of the function [K 2 e(£1, K)\ 1 in the 
complex K x plane, and their modification from Im(O) > (left panel) to the limit Im(fJ) — > 0+ 
(right panel), for < Re(f2) < Qq(K^), with Qq(K^) defined from (|A4]). The branch cuts of the 
square root and the logarithm are shown with the black dashed lines and the blue dot-dashed 
lines, respectively. The poles ±K X at the roots of e(Q,K) = are shown with the filled circles. 
The arrows show the direction of motion of the singularities in the process of Im(f2) — > 0+. The 
contour of integration over K x in (|16p is shown with the solid red line, with the arrow showing the 
direction of the contour. 

the case considered here; therefore, for any Im(f2) > the poles ±K T X are located away 
from the real axis of the complex K x plane, and thus do not lie on the integration 
contour in (1A1I) . 

The location of these singularities in the complex K x plane is shown in Figs [3H5J For 
Im(f2) > (the left panels of Figs [3H5]), none of the singularities of [K 2 e(£l,K)] 1 lie 
on, or intersect with, the real axis lm(K x ) = of the complex K x plane, along which 
the integration in flAll) is carried out; the function flAll) and, consequently, the function 
^(f2,i^||) are thus analytic functions of Q for Im(f2) > 0. However, as Im(f2) — > 0+ in the 
process of analytic continuation of K») to the real axis of complex Q plane, some of 
the singularities of [K 2 e(Vl, K)] 1 move about in the complex K x plane as shown in Figs [3]- 
[51 and may collapse onto, or cross with the real axis, thus requiring deformation of the 



integration contour in flAlj) to avoid crossing these singularities and to preserve analyticity 
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FIG. 4: (Color online) Same as Fig. but for Re(fi) > Sl (K\\) > 0. The branch cuts of the 
square root and the logarithm are shown with the black dashed lines and the blue dot-dashed 
lines, respectively. The poles ±K X at the roots of e{Q, K) = are shown with the filled circles. 
The arrows show the direction of motion of the singularities in the process of Im(r2) — > 0+. The 
contour of integration over K x in (|16j) is shown with the solid red line, with the arrow showing the 
direction of the contour. 

of C(f2, K\\). Below we consider three cases: (i) |Re(0) | < \VLq{K\\)\, (ii) |Re(0) | > \Vt {K\\)\, 
and (iii) Re(fi) = ■±.Q,q(K\\), where ±Qq(K\\) G K are found from the equation 



e(n,K)\ 



K x =0 



0, for VI, K» e 



(A4) 



For |Re(f2)| < |fio(-R"n)|, the logarithm branch cuts (1A2|) collapse onto the real axis, while 
the poles ±if£ collapse onto the imaginary axis of the complex K x plane as Im(fi) — > 0+. 
The resulting integration contour in K x is shown in Fig. [3l 

For |Re(fi)| > |J7o(iT||)|, both the logarithm branch cuts ( 1A2I) and the poles ±K X collapse 
onto the real axis of the complex K x plane as Im(fi) — > 0+. The resulting integration 
contour in K x is displaced to avoid the poles as shown in Fig. |H 

Finally, for Re(fi) = ±f2o(i^||) (the boundary between the above two cases), the two 
poles ±iC£ both collapse towards zero as Im(f2) decreases, as shown in Fig. [51 "squeezing" 
the K x integration contour between them. In the limit Im(f2) — > 0+, the poles ±K X merge 
at K x = 0, and the integration contour passes through both of them, resulting in ((Q,K\\) 
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FIG. 5: (Color online) Same as Fig. H but for Re(fi) = fio(#j|) > 0. The branch cuts of the 
square root and the logarithm are shown with the black dashed lines and the blue dot-dashed lines, 
respectively. The poles ±K X at the roots of e(fi, K) = are shown with the filled circles. The 
arrows show the direction of motion of the singularities in the process of Im(O) — > 0+. The contour 
of integration over K x in (|16j) is shown with the solid red line, with the arrow showing the direction 
of the contour; for O = ±f2o(-K"ii), this contour is "squeezed" between the poles ±K X , leading to a 
singularity in £(f2, ifn) for Im(fi) — > 0+, when the poles — K r x and K r x merge. 

being singular (non-analytic) at the point Q = ±Qq(K\\). 

The function ((Q,K\\), denned for {Q,K\\} 6 R as described above, was calculated nu- 
merically by performing integration over K x along the appropriate one of the integration 
contours shown in Figs EHU The characteristic ^-dependence of thus defined ((Q,K\\) is 
shown in Fig. [61 for a fixed value of Ku > 0. 

Appendix B: Dispersion of surface waves for Ku 3> 1 

Consider the function ((fl,K\\) defined by Eq. (|T6|) . with e(fl, K) given by (fT7|) . whose 
definition is extended to Q G 1R as described in Appendix [K\ For any given Q e R, one 
can find K» large enough so that |Q| < |O (-^||)|) where f2 (-^||) is defined by Eq. (1A4|) (this 
follows from the fact that Q{K\\) is a monotonically growing function of K\\). Hence for 
sufficiently large K«, all poles of the function l/K 2 e(Q, K) under the integral in f)16p lie on 
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FIG. 6: (Color online) The characteristic ^-dependence of ((Q,K») analytically continued to 
Im(J7) = 0, for a fixed Ku = 0.1. Re[£(fl, -Bui)] and lm[£(0, Kn)] are shown with the solid blue and 
dashed purple lines, respectively. The frequency Oo(ifii) defined by Eq. (|A4|) . at which ((£l,K\\) is 
singular, is shown with the dotted vertical line. 

the imaginary axis of the complex K x plane, symmetrically above and below K x = 0, as 
shown in Fig. [3j Therefore, the "projection" of 1/ K 2 e(Q, K) onto the integration contour 
along the real axis of complex K x plane is peaked at K x = (the point on the integration 
contour closest to the poles) and decreases away from K x = (as the distance from the 
poles increases). At K\\ ^> 1 the poles K x = ±iK\\ of 1/K 2 are located closer to the real 
axis of complex K x plane than the poles K x = ±K X of l/e(Q, K), hence the variation of the 
integrand in £(fl,K\\) with K x is mainly defined by the variation of 1/K 2 = 1/{K 2 + K 2 ) 
along the real axis of complex K x plane. Neglecting the variation of l/e(Q, K) with K x (and 
using the peak value of l/e(fl,K) at K x = 0), we can thus approximate the integral over 
K r as 



+o ° dK T r +o ° dK T 



7T 



(Bl) 



K 2 e(n,K) ~ (K 2 x + K 2 )e(^K\\) s{Sl,K()K\\ 

with e(Q, K\\) defined by (TPTj) with K — K». The corresponding approximation for ((Q, Kn) 
is then 

-l ' 



C(n,K 




i n , [n + VzKu 

1 + ^ttt r In 



Kn > 1. (B2) 



Kj 2V3Kf { \n- V3K\ 
Substituting (IB2j) into the dispersion equation (120]) for weakly damped surface waves, we 
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have 



ft , ft + v^a'i 




(B3) 



2\/3A| ft - v^i^i 
Assuming the solution of flB3p to be of the form 



ft = V3K {] [1 + 5n(K {l )], with lim |<5ft(A||)| = 0, 



(B4) 



we obtain for 8Q(K\\): 



5ft (Ay) « 2exp [2-4AJ] 



which tends to zero at large K», in agreement with the assumption flB4[) . Thus we arrive at 
the approximation (123!) for the frequency of surface waves at K»^>1. 
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